Long-Range Coulomb Effect on the Antiferromagnetism in Electron-doped Cuprates 
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Using mean-field theory, we illustrate the long-range Coulomb efi'ect on the antiferromagnetism in 
the electron-doped cuprates. Because of the Coulomb exchange effect, the magnitude of the effective 
next nearest neighbor hopping parameter increases appreciably with increasing the electron doping 
concentration, raising the frustration to the antiferromagnetic ordering. The Fermi surface evolution 
in the electron-doped cuprate Nd2-a;Cea;Cu04 and the doping dependence of the onset temperature 
of the antiferromagnetic pseudogap can be reasonably explained by the present consideration. 
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The Fermi surface (FS) evolution with electron doping 
in Nd2-xCe3;Cu04 has been observed by angle-resolved 
photoemission spectroscopy (ARPES) experiments [1,2]. 
Recently, much attention has been paid to understanding 
the physics of this phenomenon [3-8]. At low tempera- 
tures, the electron-doped cuprates are in the states of 
antiferromagnetic (AF) phase within a wide doping re- 
gion [9,10]. The FS evolution is therefore closely relevant 
to the antiferromagnetic correlation. From the frame- 
work of the Hubbard model, the doped electrons occupy 
the upper Hubbard band. On the other hand, according 
to the experimental observation, with increasing electron 
doping, the energy gap should decrease and eventually 
close up at the optimal doping, x ~ 0.14, where the AF 
phase terminates. Therefore, to interpret the experimen- 
tal results, one needs to assume the on-site interaction 
U in the Hubbard model to be dramatically decreasing 
with increasing electron doping concentration. This is a 
puzzle within the Hubbard model with constant on-site 
interaction. Some investigators have treated the doping 
dependence of U by considering some kind of screening 
[3,4,6]. ^ 

In this work, we explore the long-range Coulomb ef- 
fect (LRCE) on the AF ordering. Due to the Coulomb 
exchange effect, the LRCE results in excess electron hop- 
ping. As a whole, the excess hopping tends to frustrate 
the AF ordering. With increasing the electron doping, 
this exchange effect becomes significant, leading to the 
decreasing of the energy gap. We will see this gives a 
reasonable explanation of the FS evolution and the en- 
velope of the onset temperature of the antiferromagnetic 
pseudogap in electron-doped Nd2-a:Cea;Cu04. 

We start with the following two-dimensional square 
lattice model 



H 



(1) 



where tij denotes the hopping energy of an electron be- 
tween the lattice sites i and j, {cia) represents the 
electron creation (annihilation) operator of spin-CT(= ±1 



for up and down spins, respectively) at site i, ni„ — 
cl^Cia, and rii = riii + "■U electron density oper- 

ator. Besides the on-site Hubbard repulsion U, we take 
into account the long-range Coulomb interaction in the 
third term. Such a similar type of Hamiltonian has been 
used to investigate the LRCE on the d-wave pairing for 
the hole-doped case [11]. 

In the hopping term, besides the nearest neighbor (n.n) 
hopping, the other hopping processes within a range need 
to be included. The essential role of the next n.n hop- 
ping for the validity of the single-band Hubbard model 
for describing the cuprates has been investigated by com- 
paring it with the two-band [12] and three-band Hubbard 
models [4]. The particle- hole asymmetry in the cuprates 
can be understood by taking into account the next n.n 
hopping [13,14]. By including the next and third n.n hop- 
ping parameters in the types oft — J and Hubbard mod- 
els, numerous studies have been carried out to explore 
the properties of the electron-doped cuprates [3-8,15-18]. 
We here take two more additional parameters for the 
fourth and fifth hopping, each of which is smaller than 
the formers. In terms of the n.n hopping parameter, 
i(i 0) = t, the values of other 4 parameters are given 
as = -0.325i, i(2,o) = 0.17i, i(2,i) = -O.Ult, and 

^(2,2) — — 0.07t, where the subscripts instead of ij denote 
the components of a vector of length |i— j|. The values of 
the next and third n.n hopping parameters are approx- 
imately the same as that in the literatures [3-8,15-18]. 
From the hopping term, the single-particle dispersion e° 
is given by 

Cfc = -'2t{Cx + Cy) - U(^i i)CxCy - 2<(2,o)(c2a; + C2y) 
-4t(2,l)(c2xCy -I- C2yCx) - ^t(2,2)C2xC2y 

with cix — cos{lkx). 

For the long-range Coulomb interaction Vij , we take 
the following form 



yiexp(-|i-j|/d)/|i-j| 



(2) 



where the length scale is in unit of the lattice constant. 
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Here we set Vi = l.Ot and d = 4 for the present calcula- 
tion. For the on-site U, we use U = A.St, which is within 
the range of the doping-dependent interaction adopted in 
the existing calculations [3,4,6]. The strength of Vi with 
Vi/U w 0.21 is a typical one [11]. 

By the mean-field (MF) approximation, from the long- 
range Coulomb term, we get the exchange part of the 
self-energy as 

^^ = -^E<|k-k'l)(4'TC'^'T) (3) 

k' 

where TV is the number of lattice sites, ti(|k — k'|) is the 
Fourier transform of , and (■ ■ ■) denotes an average. 
This exchange self-energy is equivalent to a hopping en- 
ergy in real space. The corresponding hopping integral 
is given as if^ = ''^iji'^i'^n)- ^^'^ significance of this 
exchange effect and its doping dependence will be dis- 
cussed later. With the contribution from the exchange 
self-energy, the effective single-particle dispersion is given 
hyek = el + ^l. 

At low temperatures, the strong on-site Coulomb in- 
teraction leads to the AF ordering. In the case of the 
electron doping under consideration, the AF ordering is 
a commensurate spin-density wave. The order parameter 
is given by A = U{ni{Q) — n^{Q)) /2N , where n„{Q) is 
the Fourier transform of the electron density of spin-ci at 
wave vector Q = {tt,tt). By the MF approximation, the 
energy spectrum of the quasiparticle is then given by 

= (efc + ek+Q ± Ek)/2, (4) 

where the -I- (-) sign refers to the upper (lower) Hubbard 
band, and ~ {ek ~ e/c+Q)^ + 4A2. The parameter 
A and the chemical potential /x are self-consistently de- 
termined by 

k 

k 

where / is the Fermi distribution function, and n is the 
electron density. In terms of doping concentration x,n — 
\ + x. 

To analyze the ARPES observations on the FS evolu- 
tion, one usually considers the spectral density occupied 
by the electrons. By the MF theory, this spectral density 
is given by 

A<(fc, = j{uj){v?(k)b(uj - El) + v^(k)b{uj ~ El)\ 

(6) 

with u^{k) - 1 - = [1 + (efc - efc+Q)/£;fc]/2. Fig. 

1 shows the Fermi-level electron distributions in momen- 
tum space obtained by integrating A^ik^iS) in a energy 



window \io — [l\ < 0.15t at temperature T = O.lt. At 
low doping, the FS (in the first quadrant of the Brillouin 
zone) appears as two small pockets centered at (tt, 0) and 
(0,7r). With increasing doping, these pockets extend to 
larger squares. Such an evolution of the pockets reflects 
the occupation of the doped electrons at the upper Hub- 
bard band. On the other hand, a small pocket begins to 
form around (7r/2,7r/2) at x ~ 0.1 and grows up with 
increasing doping. The formation of this pocket stems 
from the contribution of the lower Hubbard band; with 
increasing doping, the energy gap decreases and mean- 
while the lower Hubbard band shifts up toward to the 
FS. At high doping, the energy gap closes up, therefore 
the FS is a single curve centered at (tt, tt) {x = 0.18). 
Clearly, the present calculation reproduces the results of 
the existing calculations based on the doping-dependent 
U Hubbard model [3,6] and is in agreement with the ex- 
perimental observations [1]. By the t — t' — t" — U model, 
a constant U actually results in a nearly constant en- 
ergy gap in a wide region of the doping concentration, 
and hence cannot explain the FS evolution. Though U 
is constant in the present calculation, but because of the 
LRCE, the energy gap decreases remarkably with increas- 
ing doping. We will discuss this problem later again. 



(0,0) X = 0.04 ^ 


x = 0.l 


i 

K 

jc = 0.14 


:c = 0.18 



FIG. 1. (Color online) Maps of Fermi surface at various 
doping concentrations at temperature T — O.lt. Highs are 
denoted by red and lows by blue. 

In Fig. 2, we depict the MF Neel temperature 
(solid line with circles) as a function of the doping con- 
centration X, and compare it with the result of the ex- 
tended Hubbard model (EHM, in which the LRCE is not 
included) and the experimental data. It is seen that the 
difference between the present calculation with LRCE 
and the EHM is remarkable. By the EHM, the depen- 
dence of Tjv on X is very weak in the region < a; < 0.2. 
In contrast to the EHM, Tjv given by the present cal- 
culation with LRCE decreases distinctly with increasing 
doping. In particular, it drops sharply at x ~ 0.141. 
This doping concentration corresponds to the quantum 
critical point at which the zero-temperature AF transi- 
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tion terminates and is consistent with the experimental 
observation [9,10]. Also, the envelope of with LRCE 
is in fairly good agreement with the experimental results 
for the pseudogap Apg and its onset temperature T* [10]. 
The MF Neel transition is associated with the appearance 
of the local magnetization, and hence the comparison be- 
tween the MF Tn and T* or Apg is meaningful [4] . The 
MF order parameter is a measure of the local order of 
the real system with long-wavelength fluctuation, similar 
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FIG. 2. (Color online) Normalized MF Neel temperature 
Tjsi{x)/Ti^{0) as function of the doping concentration x. The 
solid line with circles is the present calculation taking into 
account the LRCE. The dashed line with squares is the re- 
sult of the extended Hubbard model without the long-range 
Coulomb interaction. The pseudogap Apg and its onset tem- 
perature T* (both of them normalized by Ao = 4:60K) are 
the experimental results [10]. 

Why does the LRCE lead to Tjv decreasing with in- 
creasing doping? To answer this question, we here 
consider the excess hopping parameter tfj due to the 
Coulomb exchange. This quantity can be written as 



,x 



2N 



^(cj,|Cfe|)5/(/c) 



(7) 



where Si{k) ~ cos{lxkx) cos{lyky) + cos{lyk,j;) cos{lxky) 
and again / = {lx,ly) represents a vector of length |i — j|. 
The doping dependence of originates from the Fermi 

distribution (c[.|Cfc|). At low temperature, the contribu- 
tion to the integral in Eq. (7) comes from the Fermi area. 
The Fermi surface varies with changing doping concen- 
tration. Especially, with increasing doping, the Fermi 
area increases considerably in the regions near the points 
(tt, 0) and (0, tt) since at which the upper Hubbard band 
reaches its minimum and where the energy dispersion is 
nearly flat. Therefore, the variation of tf stems predom- 
inantly from the integral in these regions. On the other 
hand, with the change of the Fermi surface, the factor 



Si{k) varies slowly for small I, but rapidly for large 
It is then easy to see the variation of tf with increasing 
doping at low temperature: for example, i^-^ is almost 
unchanged for 5(i.o)(7r, 0) = 0; tf^i has a large vari- 
ation because of S'(i^i) (tt, 0) = —2. Since the original 
hopping parameter is negative, the magnitude of 

the effective hopping parameter + t^i 2) should be 
enhanced with increasing doping. For large I, because 
of the cancellation from the destructive factor Si{k), not 
only the variation, but also the magnitude of tf are neg- 
ligible small. At high temperature, since the integral in 
Eq. (7) is taken over a spreading area wider than the 
Fermi area, the behavior of tf is not so intuitive. Fig. 
3 exhibits the variations of the magnitudes of the effec- 
tive hopping parameters ti + tf as functions of the doping 
concentration x at temperature T — 0.53i. The quantity 
Sti{x) is defined as 



6tiix)^[tf{x)-tfmsgn[ti+tfm- 



(8) 



From Fig. 3, it is seen that except for (5t(i_o) that is nearly 
constant, the other parameters vary approximately lin- 
early with X. Especially, among those effective hopping 
parameters only the magnitude of + t^^ = t' in- 
creases appreciably with increasing x. This is consistent 
with the above analysis. It is know that a large ratio t' /t 
can destroy the AF instability at weak U [21]. This is 
also true for strong on-site interactions. For large t' /U , 
the next n.n. AF coupling constant is J' — At''^ /U, which 
leads to the frustration of AF ordering in the square lat- 
tice. It is therefore clear that the AF order parameter 
and the transition temperature decrease with increasing 
doninff concentration. 
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FIG. 3. (Color online) Variation of the effective hopping 
parameter 5ti as function of the doping concentration x at 
temperature T = 0.53f. A pair of numbers with brackets 
denotes the vector I = {lx,ly)- 

Shown in Fig. 4 is the parameter tf as a function of 
distance at the MF = 0-39t at a: = 0.1. It is seen that 
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only the magnitudes of the first 4 parameters are appre- 
ciable; from the fifth to the one at distance 4, all of them 
are very small; after that all other parameters are negli- 
gible small. In a wide region on the temperature-doping 
concentration phase diagram, the parameter tf behaves 
almost the same as shown in Fig. 4, with only a visual 
change in the next n.n hopping parameter. From the be- 
havior of the parameter tf , the range of the predominant 
Coulomb exchange effect seems to be shorter than 3. Of 
course, these results are obtained for the particular sets 
of the input parameters (t('s, U, V\ and d). But these 
parameters are reasonable. In the present calculation, all 
the innut narameters were not fine-tuned. 
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FIG. 4. The hopping parameter due to Coulomb ox- 
change as function of the distance i at T = 0.39f and x = 0.1. 

In summary, we have investigated the long-range 
Coulomb effect on the antiferromagnetism in the 
electron-doped cuprates using the mean-field theory. Due 
to the Coulomb exchange, the magnitude of the effective 
next nearest neighbor hopping parameter in especial in- 
creases appreciably with increasing the electron doping 
concentration. This leads to stronger frustration to the 
AF ordering at higher doping concentration. Therefore 
the transition temperature decreases with increasing dop- 
ing. Consequently, the AF phase terminates at a doping 
concentration x « 0.14 in consistent with experiments. 
The present calculation gives a reasonable explanation 
to the doping dependence of the onset temperature of 
the AF pseudogap as well as to the FS evolution in the 
electron-doped cuprate Nd2-xCcj;Cu04. 
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